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SUMMARY 


fjv 


In 1879 G. W. Hill, using his differential equation system, found a 
series of direct lunar orbits, which represent a simplification of the 
restricted three-body problem (neglecting solar parallax and eccen¬ 
tricity). If T is the period of revolution time, and T = 2mn is fixed, the 
initial conditions, <f 0 (m), rj 0 (m) , and (f 0 = h 0 = °)» are dependent upon m 
and give periodic solutions. Hill obtained these functions in the form 
of power series, with respect to m, which converges rapidly for small 
distances of the satellite; and he gave their first coefficients. 


In the present discussion, Hill's coefficients are successively de¬ 
termined from linear equations and are in the form of closed trigono¬ 
metric expressions. Although this method—which does not extend 
beyond the third order—is only an interesting attempt at achieving the 
same goal via a different mathematical route, it does not extend as far 
as Hill's method (which may be carried out to the fifth or sixth order 
without difficulty). 
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REMARKS ON HILL’S LUNAR THEORY. 
PART II 


by 

Karl Stumpff* 

Goddard Space Flight Center 


INTRODUCTION 

In his famous writing of 1878, G. W. Hill^ found a series of direct lunar orbits as solutions of 
his differential equation system. This series consitutes a simplification of the probleme restreint 
(by discounting the solar parallax). Hill found that if T is the (synodic) circuit time of the satellite 
and T = 277 m is fixed, the initial conditions are £ 0 (m), 7? 0 (m), (<f 0 = tj 0 = o), and are dependent upon m 
as a parameter and give such periodic solutions. He obtained these functions in a power series with 
respect to m, which converge rapidly for small distances of the satellite (short circuit times, small 
m), and gave their first coefficients. In the following examination, another method will be presented 
which leads to the same goal, but which permits these coefficients to appear as finite trigonometric 
expressions. 


HILL’S DIFFERENTIAL EQUATIONS 


Hill’s differential equations 
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where r 2 = <f 2 +7? 2 , are based upon the rectangular coordinates 
of a massless satellite in reference to its central body (planet) 
and in a coordinate system which revolves at the same angular 

S — 


o 

CL 


velocity, n = 1, as the planet which orbits the sun (Figure 1), The 

, , , a . Figure 1 —Geometric representation 

location of the sun S is therefore assumed to be in thedirection of D f a mas5 | e$ s satellite sofa planet 

the negative f-axis, at any desired distance, so that the parallax P. 

•NAS-NASA Research Associate; Professor Emeritus, Gottingen University. 

tHill, G. W., "Researches in the Lunar Theory,” in The Collected Works of George William Hill, Vol. !, Washington, D. C.: Carnegie 
Institution, 1905* 
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of the sun can be discounted. Under these circumstances, Equations 1 apply for the satellite S. 
Periodic solutions can be expected if we let the satellite start from a point (<f 0 , 0) on the f-axis with 
a velocity (0, t? 0 ), whereby the initial (positive) velocity must occur perpendicular to the <f-axis and 
t] 0 will depend to a certain extent upon f 0 , Every periodic solution of this sort is symmetric to both 
axes, since Equations 1 do not change when iis exchanged with and 77 with -77. If T = 277 m is the 
circuit time in such a periodic orbit, then the total of these orbits also can be characterized by de¬ 
termining <f 0 = <f 0 (m), i) 0 = 7? 0 (m) as functions of the parameter m. 

If r = 0 is the moment at which the satellite crosses the positive £-axis at right angles, then the 
periodic motion has the form 

£(t) = Aj cos vr + A 3 cos 3vt + A $ cos 5 vr + — , 

77(7-) = Bj sin vt + B 3 sin 3vr + B s sin 5vr + • - * , 

where v = i/m = 2jt/t 0 For circuit times as brief as desired (m -0) —hence for very close orbits —the 
perturbations by the sun become as small as desired. The above equations (Equation 2) now change 
into 

£(r) = a cos vr, rjir) - a sin vr , (3) 

or into the form of uniform circular motion, for circles are the only undisturbed orbits which are 
symmetric to both axes. If we differentiate Equation 3 any number of times, we find, for this bound¬ 
ary case and for r = 0, 



„2 n ^(2n) = (-l) n f 0i 


n 2«-1^2n-l = . 


The above equations mean that, for m - 0, the quantities 


^0 


1 . 01 = 


7} ° 

+ ffl “—-* 1 , 


+ m 4 


^0 ( 


( 4 ) 


1 , 0. = 


, V <> 

-m 3 -J— - 1 , 
s 0 


a 6 ” _m6 


#n (6) 
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all approach unity. 
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If for the sake of abbreviation we set l/r 3 - ^ then the differential equations (Equation 1) can 
be written: 

} = 2y - £(/x-3) ; and V? = “2 g - w . 

If we differentiate these equations any number of times and if we consider that uneven derivatives of 
<f and Mand even derivatives of v disappear for r = 0, we then obtain the system: 


- m 

+ m 4 

~ m 6 

^0 


^0 

+ m 5 
^0 

-m 7 

If these equations are 
become: 


fo 27? o ^0(^0 1 

^o <4> = 2 *>o ■ |/o^o~ 3 ) + fo^o] ’ 

^ o <«) = 2r, 0 (S> - |> 0 < 4 >(m 0 -3) + 6^ 0 Mo + ^o^d (4) ] 

% = ~ 2^o “ ^0 M - 

^o <5; = " ^0 4) " (**>> + 3 ^) • 

^(7) = -2^ 0 («) - (»,,(*> M + 10 v 0 Jj. + St) 0 m <4> ) . 

multiplied by the factors to their left and if Equation 4 is substituted they 


“ 2m/3 1 

“ 3m 2 + 4> , 

Pi 

= * , 


-2^3 

- 3m 2 a 2 + 0v 2 , 


= “ 2ma 2 

+ <^i/j 

-2m/3 s 

- 3m 2 a 4 + , 


= - 2ma 4 

+ ^.v 3 

-2nvfl ? 

~ 3m 2 + q6v 6 , 

Pi 

= -2ma 6 

+ <£v 5 
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where 


4> - m 2 fj.Q and \jt = p i 




m f ; 

g o 


and, with 
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V 2 = a 2 + M 2 


"N 


v 3 - /3 3 + 3/2,^ 


V 4 " a 4 + 6 ^2 tt 2 + ^4 


( 8 ) 


v 5 = p s + 10m 2 P 3 + 5^„ /3j 


V 6 " a « + 1S ^2 “4 + 1S ^4 a 2 + ^6 • 


If the quantities m 2 > ^ 4 * ' are eliminated from Equation 8 by forming the even derivatives of m and 

substituting their values for r = 0 into Equation 7, the following system of recurrence formulas is 
obtained for v. ; 

*1 = /»! ■ 

V 3 = ^3 " 3a 2 V l + 3 - fl l V 2 - 


£5 " Sa 4 V l + 10 ^3 V 2 ~ 10a 2 ^3 + ^4 


M9) 


-(1 + 1) a 2 + (1 +2)/3,v 1 , 

-(1 + 1) a 4 + {3 + 4) Vj - (3 + 6) a 2 v 2 + (1 + 4) fi t , 

-(l + l)a 6 + {$ + 6)/3 s v t - (10 + 15)a 4 z/ 2 + (10 + 20) i> 3 - {5 + 15)a 2 z/ 4 + (1+6 )P x v $ , 


In these formulas we again recognize the binomial coefficients of uneven order in the coefficients 
of v { with uneven indices, while the coefficients of v i with even indices are made up of the binomial 
coefficients of uneven and subsequent even orders. In all cases, the signs are alternating and the 
sum of the coefficients is unity. However, this means that all v., just as the a. and /3., approach 
unity for m-0. The same also applies for 4> and 0, as can be seen from Equation 5. With the help 
of Equations 5 and 9, a. and fi. can be determined one after another, in the order of increasing indi¬ 
ces (p l9 a 2 , p 9 a 4 , • • *), as functions of m, 4> y \p. 
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The orbit of a satellite which originates from (f 0 , 0) at time r - 0 with the velocity (0, t ? 0 ) will 
be periodic if it reaches the point (-f 0 , 0) at r = (1/2 )T = 77 m. Therefore the periodicity conditions can 
be written: 

£(77111) = - £ 0 , 17(77111) = 0 , ( 10 ) 


or, if we set up the Taylor series, 


£ (77TT1) 


-fo = fo 


(77m} 
2 ! 


(77m ) 4 
4! 


+ 


17(77111) - 0 


77 m 77 0 


( 77TT1 ) 

~3T 


^0 


(77m } 5 

“TT“ 


’V 


( 5 ) + 


If we divide by £ 0 and substitute the quantites of Equation 4, we can then write 

7T 2 77 4 77 6 

f(m) - 0 - 2 - a 2 — + a 4 -gr + * • • 

g(m) = 0 = + /3 s 3T - ••• • 




^ (ID 




On the other hand, 


f (m) 


g(m) 


“ go 


f o m + 


go m 


» 17,2 
f 0 TT 


m 

2 ! 


+ f 


0 3! 


m jnl 
go 3! 


* 0 


- 0 


y ( 12 ) 


if the index 0 is based upon m = 0 and the derivatives according tom indicated by primes. The Equa¬ 
tions 12 are identities, since they are fulfilled for any given value of m. Thus the requirement of 
periodic orbits is given by 


f 


o 


f 0 




0 , 


g 0 go go 


0 . 


(13) 


In general, then, if n is any positive whole number or zero and if a/"), are understood to be the 
n th derivatives of the quantities a. , p i with respect to m for m = 0 , then 
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2 ~ 2 ! + TT ” "6T + * ‘ ‘ - 1 + cos it - 0 , 


ITT + “ST ~TT + * * ’ = sinyr - 0, 


(14) 


since all a., p i assume the value 1 for m - 0, and furthermore 


f 0 (B) = 0 


a i Cn) TT + a 4 (n> IT" a 6 (n) TT 


g 0 (n> = o = -/8j(«>-jr +/3 S ( n > - 


(15) 


The Equations 14 are identically fulfilled and only show that the arrangements were correct. 
For additional computations we shall make an analysis of Equations 15 for n = 1 explicitly; i.e., 


V 


= o = - 


a i it + a A ~ a 


77 

« 61 


77^ 77^ 

0 = fi t ' TT ‘ ^ 3f + TT 


(16) 




We have shown above that it is possible to represent the quantities a i> /?. as functions of m, <p and 
Since <f> and \p assume the value 1 for m = 0, and if we set 


<£(m) = 1 + ajm + a 2 "JT + a 


3 3! 


(17) 


^ (m) 


1 + b.m + b, -~r + b 


2 2! T °3 3! 


it follows, for m = 0, that 


*0 (n) - a n - ^0 (n) - b n ■ (18) 

From Equations 5 and 9 we obtain, one after another, the equations: 

^ /V = r . 

a .2 ~ ~ 2/3j _ 2nv3j' - 6m + <$>' , 
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P 3 ' = - 2 a, - 2 ma 2 ' + <£', 3 , + 0/3,' , 

a 4 ' = - 2/3 3 - 2m/3 3 ' - 6ma 2 - 3m 2 a 2 ' + + <*>(-2a 2 ' + 6/3,£,') , etc. 


If, in these and the following formulas, we set m = 0 and consider that the a., /3., v., <t>, and 4> go 
towards 1 and the derivations 4>', <P‘ towards a,, b,, then, after some brief calculations, we get 


a 2 3 l 2 ’ 

“«' - 6b r ( 3 1 ~ 2 ) ' 

- 48b l " 21 (a, -2) , 
a,' - 234b, - 113 (a, - 2 ) 


Pi - b . - 

N - b i + ( a , - 2) • 

Pi - 19b, - 7(a,-2) , 
Pj - 109b, - 51 (a, -2} 


If we substitute these expressions into Equation 16, we get two linear equations: 

Abj - B( ai -2) = 0 , 

Cb x - D(a 1 - 2) = 0 , 


(19) 


the coefficients of which are represented by the 77-series: 


C = 


D = 


= 64 

r ■ ■ 

77 2 

77 4 

" 2! 

+ TT 

= 77 - 

77 3 

3! 

77 3 


" 3! 

+ 7 “ 


48 TV + 234 


8T 


7T J 77 

+ - 109 -JT + 


5! M 7! 


7 



This series can be easily summed, for from 


77* 77 7T 9 77° 


1 2! + 4! ITT + KT 


1 _ 1 77 2 77 4 7T W 77 

"7 cos 277 - “j - 2 "JT + 8 "^T “ 32 “gT + 128 


7® 

6T 


8 

ST 


_ 77 2 77 4 77 ® 77 ® 

277sm77 - 4 "2T ” 8 ~^T + 12 "gr ~ 16 "gT + 


we obtain by summation: 


cos 77 + -y cos 277 + 277 sin 77 


3 ( 77 2 77 4 77® 77 ® \ 

T + (tT + 4! - 21 gr ♦ 113 3 T - 


3 

= -j + B 3 " 2 


Similarly we represent A, C,D by closed expressions: specifically, 


3 + 2 cos 77 + cos 277 + 377 sin 77 - -4 


_ 3 1 

B - “ 2 + COS 77 + ycos 277 + 277 S in 77 - ”2, 


C - 2 s in 77 + s in 277 - 3?7 cos 77 - 377 } 


s in 77 + ~2 s in 277 ~ 277 cos 77 - 277 


r ( 20 ) 


Since the determinant of the homogeneous system (Equation 19), BC-AD = 27?, differs from zero, it 
follows that 


a, - 2, b, - 0 . 


( 21 ) 


After a t , bj are known, we can determine a 2 , b 2 in a corresponding manner from the pair of 
equations f" = 0, g" = 0. From these a pair of linear equations results: 


Ab 2 - Ba 2 = X 2 
Cb 2 - Da 2 = Y 2 


(22) 


8 



in which A, B, C, D are defined as in Equation 20, and in which X 2 , Y 2 represent new ^-series. The 
coefficients of this v -series can be determined numerically by using Equation 21. With the values 
(Equation 20) forA,""D, we find the general expressions 


Ab - Ba 

n n 



Cb - Da n 

n n 



> 


(23) 


where X n , Y n are known, if all a i , up to 1 = n - 1 are present numerically. 
By carrying out this relatively easy analysis to the order n = 3, we obtain: 


Xj = 3 “ 2 cos it - cos 2 t7 - 477 sin it - ~ 2B - +4 


19 3 3 

X 2 = - 3 + COS 77 + COS 277 - -g- COS 377 “ ~ 2 ~ Tf Sin 77 - “ 4 


X 3 - — 2 cos 77 —2 cos 277 - 3 cos 3 t 7 - 3677 sin77 - +10 


Y, - - 2 s in 77 - s in 277 + 4 tt cos 77 - - 2 D - - 477 , 


19 3 3 3 

-q* sin 77 + sin 277 -^"sin37T + -y cos 77 - 2 77 


Y 3 - -22sin77 — 5 "sin 277 - 3sin377 + 3677 cos v - ~ 36 tt 


(24) 


The first coefficients of the development (Equation 17) are therefore known, and we obtain 


<£(m) - 1 + 2m + 9 2 1 + 57 3 ! + ’ ' * 


and 


11 m 2 m 3 

\ p [ m ) - 1 + 2 2* + ^6 3? + 


But according to Equation 6, we have 


irr 

<5 3 ’ 
* o 


^0 

m 

S n 
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so that we can also calculate 4> and ^ from Hill's power series for £ 0 and 7? 0 . This calculation yields: 


<£( m) 




m 2 m 3 1649 ra 4 1235611 m 5 

1 + 2 m + 9 -jr + 57 -jr + *T + 324 TT + '" . 


11 m 2 m 3 3857 m 4 23153 m 5 

1 + TT + 26 TT + IT" TT + — IF” T T + 


and completely verifies our analysis. 


(25a) 

(25b) 


CONCLUDING REMARKS 

The foregoing method of obtaining Hill’s coefficients is not intended to be more than an interest¬ 
ing attempt to achieve the same goal via a completely different mathematical route. This attempt is 
interesting for two reasons: (1) Because the coefficients of the series (Equations 25a and 25b) are 
determined successively from linear equations (Equation 23); and (2) Because they are in the form of 
closed trigonometric expressions. It must not be overlooked, however, that this method — elegant 
though it may be — does not lead as far as Hill’s, the developments of which can be carried out to the 
fifth or sixth order without difficulty. Here we could hardly pass beyond the third order — at least 
not with the method used, which is, however, probably not the only possible one. The reason for this 
is as follows: 

1. The method of determining the coefficients of the trigonometric formulas (Equation 24) may 
be demonstrated by using the example of X 2 , for which we obtain the n -series 


X 


2 





77 ° 

2190 gT + 


(26) 


Then the arrangement 


= a + b cos 77 + c cos 277 + d cos 3?7 + e77 sin 77 


- a + b 





( TT 2 77 * \ 

+ c ! i - 4 -jr + 16 ~4f - ■ ■ • J 

/ 77 2 77 * \ 

+ d fl - 9 ~ 2 "T + 81 “ 4 T ” 


+ 


e 







> (27) 
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leads to the equations: 


a + b + c + d 

= 0 , 

- b - 4c - 9d + 2e 

= -6 , 

b + 16c + 81d - 4e 

= ‘6 , 

-b - 64c - 729d + 6e 

= 198 , 

b + 256c + 6561d - 8e 

= - 2190 




r (28) 


J 


the number of which can be expanded as far as the coefficients of the v-series are known. Only if 
the first five of these equations yield the coefficients of the trigonometric arrangement (Equation 27), 
will this formula solve the problem exactly; however, when the coefficients — a, • • -e — calculated by 
means of it also fulfill the other equations (Equation 28), a demand which could only be satisfied if it 
were possible to find a law (and then to prove it) by which the terms on the right side of Equation 28 
— the coefficients of Equation 26 — would increase. 

In the case of the developments (Equation 24), it was possible to check the validity of the form¬ 
ula with at least one of the supernumerary equations and to thereby prove it. Furthermore, in this 
case, it is not difficult to find the law of the progression of the coefficients. The agreement with the 
results of Hill also vouch for their correctness. 

2. The difficulties of extending this method begin with the determination of X 4 , Y 4 , for the nu¬ 
merical values of the coefficients of the7T-series (Equation 26) increase extraordinarily, and be¬ 
cause until now it has remained uncertain which formula (Equation 27) would achieve the right goal. 

A trial of the formula 


X 4 = a + b cos 7 t + cos 2rr + d cos 377 + e cos 477 + f77 sin tt 

and a corresponding one for y 4 has not succeeded, probably because the formula requires still more 
terms — at least one with cos 5n, 

RECOMMENDATIONS 

Further experiments with this mathematical problem could be useful and interesting, and we 
shall conclude this report with some suggestions: 

1, The expressions of Equation 24 can be given in a simpler form. For instance, we could 
write 


Xj = 2(1-costt) + 2sin 2 77 - 47rsin77 t 
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and by the substitution of the "c-functions ,T known from the two-body motion,* 


c 0 = cos A. , c l - (sin\}/\ t c 2 = (l-cos\)/\ 2 , c 3 = [k~ sin\)/\ 3 , , 

we then obtain, for \ = n: 

X 1 = wa [ 2 ( C 2 + C l 2 )" 4c l] ’ 

x 2 - ^[-2(c J+ c*) +fc 1 (c # c 1 -l)] , 

X 3 = 7T* + Cj 2 ) + 12 Cj (c 0 cj-3)] , 

and 



2. Noticeable and displeasing is the fact that the expressions of Equation 24 contain a rather 
large number of disappearing terms, which are necessarily created by the algorithm which leads to 
them. It should be possible to modify the process in such a way that these superfluous and compli¬ 
cating terms can be avoided at the beginning. Here is a mathematical problem that is not only chal¬ 
lenging, but its solution would be very useful in celestial mechanics wherever the problem of seeking 
initial conditions for periodic orbits exists. 


•Stumpff, K., "Calculation of Ephemerides from Initial Values," NASA Technical Note D-1415, 1962. 
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